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Abstract. The first LHC data on high transverse momentum hadron and dijet spectra in PbPb collisions at 
center-of-mass energy 2.76 TeV per nucleon pair are analyzed in the frameworks of PYQUEN jet quenching 
model. The presented studies for the nuclear modification factor of high-pr hadrons and the imbalance in 
dijet transverse energy support the supposition that the intensive wide-angular ( "out-of-cone" ) medium- 
induced partonic energy loss is seen in central PbPb collisions at the LHC. 
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1 Introduction 



The production of high transverse momentum hadrons 
and jets is one of the important probes of hot quark-gluon 
matter (QGM) created in ultrarelativistic heavy ion colli- 
sions. The energy loss of energetic partons, so called "jet 
quenching" , is predicted to be very different in cold nu- 
clear matter and in hot QGM, and leads to a number of 
phenomena which are already seen in gold-gold collisions 
at RHIC (v^nn = 200 GeV), such as suppression of high- 
er hadron spectrum, medium-modified angular correla- 
tions, azimuthal anisotropy, etc. (see, e.g., recent reviews 
[1, 2, 3,4,5 and references therein). The current lead- lead 
collision energy at LHC, \/%N = ^.76 TeV, is a factor 
of ~ 14 larger than that in RHIC, thereby allows one 
to probe new frontiers of super-high temperature and (al- 
most) net-baryon free Quantum Chromodynamics (QCD). 
It is expected that at such ultra-high energies the role of 
hard and semi-hard particle production may be significant 
even for the bulk properties of the created matter. 

The first striking results of LHC heavy ion program 
are the momentum-dependent suppression of charge high- 
Pr particle production [5] and the large imbalance in dijet 
transverse energy [7,8 . It was supposed in [9] with simple 
kinematic arguments that a frequency collimation mech- 
anism in medium-modified parton showering can account 
for these observables. The calculation of dijet asymmetry 
due to multiple scatterings of shower partons to large an- 
gles has been done in [TO] . The dijet asymmetry was also 
computed recently with MARTINI model in [TT] . The de- 
coherence of parton showering in a dense medium was con- 
sidered in |12j . Some numerical predictions for jet broad- 
ening in PbPb collisions at the LHC can be found in [13] . 
In present paper, the jet quenching pattern at the LHC is 
analyzed in the frameworks of PYQUEN partonic energy 
loss model [2] , and the related interpretation of the data 
is discussed. 



2 PYQUEN partonic energy loss model 

A number of Monte-Carlo models to combine a pertur- 
bative final state parton shower with QCD-medium ef- 
fects and simulate the jet quenching in various approaches 
have been developed in recent years [T31IT51ITB1I 17U18IIT9] . 
Our model PYQUEN (PYthia QUENched) QI] was con- 
structed as a modification of the jet event obtained with 
the generator of hadron-hadron interactions PYTHIA_6.4 [20] . 
The details of PYQUEN physics model and simulation 
procedure can be found in |14j . the main features of the 
model being listed only very briefly below. 

The approach to the description of medium-induced 
multiple scattering of hard partons in PYQUEN is based 
on the accumulated energy loss via the gluon radiation be- 
ing associated with each parton scattering in the expand- 
ing quark-gluon fluid and includes the interference effect 
(for the emission of gluons with a finite formation time) 
using the modified radiation spectrum dE/dx as a func- 
tion of decreasing temperature T. The model takes into 
account the radiative and collisional energy loss of hard 
partons in longitudinally expanding quark-gluon fluid, as 
well as a realistic nuclear geometry. The radiative energy 
loss is treated in the frameworks of BDMS model [2X1122] 
with the simple generalization to a massive quark case 
using the "dead-cone" approximation [33] . The collisional 
energy loss due to elastic scatterings is calculated in the 
high-momentum transfer limit [24,25,26 . The strength of 
the energy loss in PYQUEN is determined mainly by the 
initial maximal temperature T™ ayi of hot matter in central 
PbPb collisions (i.e. an initial temperature in the center 
of nuclear overlapping area at mid-rapidity). The energy 
loss depends also on the proper time To of QGM forma- 
tion and the number Nf of active flavors in the medium. 
The parameter values To = 0.1 fm/c and Nt = (gluon- 
dominated plasma) are used for the simulations presented. 

Another important ingredient of the model is the an- 
gular spectrum of medium-induced radiation. Since the 
detailed treatment of this is rather sophisticated, the sim- 
ple parametrizations of the gluon distribution over the 
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emission angle 9 are used. There are two main options in 
PYQUEN. The first one is the "small-angular" radiation, 



dN 9 
~d9 



cx sin t) exp 
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(1) 



where 6>o ~ 5° is the typical angle of the coherent gluon 
radiation as estimated in [27] . The second one is the "wide- 
angular" radiation, 



dN9 
~~d9~ 



oc 1/9 , 



which is similar to the angular spectrum of parton show- 
ering in a vacuum without coherent effects [35]. In this 
paper, in order to illustrate the sensitivity of experimental 
observables to various jet quenching scenarios, the third 
(rather extreme) option "extra wide-angular" radiation, 



dN3 
~d9~ 



i/Ve 



(3) 



is also considered. One can assume that the physical mean- 
ing of (extra) wide-angular radiation could be the presence 
of intensive secondary rescatterings of in-medium emit- 
ted gluons [in]. The collisional energy loss is considered 
in PYQUEN as "out-of-cone" loss (i.e. an energy is "ab- 
sorbed" by the medium, because the major part of "ther- 
mal" particles knocked out of the dense matter by elastic 
rescatterings fly outside a typical jet cone [27]). 

The event-by-event simulation procedure in PYQUEN 
includes the following steps. 

— The generation of initial parton spectra with PYTHIA 
and production vertexes at the given impact parame- 
ter. 

— The rescattering-by-rescattering simulation of the par- 
ton path in a dense zone and its radiative and colli- 
sional energy loss. 

— The final hadronization according to the Lund string 
model for hard partons and in-medium emitted gluons. 

In presented studies, the effect of nuclear shadowing 
on parton distribution functions, calculating the nuclear 
modification factor Raa for charged particles at px > 
10 GeV, is taken into account. The procedure is simi- 
lar to one applied in HYDJET and HYDJET++ event 
generators [29j . It uses the impact parameter dependent 
parametrization obtained in the frameworks of Glauber- 
Gribov theory [30] . The nuclear shadowing correction fac- 
tor S{ri,T2,x\,X2,Q 2 ) depends on the kinematic variables 
of incoming hard partons (the momentum fractions 
of the initial partons from the incoming nuclei and the 
momentum scale Q 2 of the hard scattering), and on the 
transverse coordinates 7*1,2 of the partons in their respec- 
tive nuclei. This initial state effect reduces the number of 
partons in the incoming hadronic wave-function of both 
the nuclei and thus reduces the yield of high-p^ hadrons 
from jet production (in addition to the medium-induced 
partonic energy loss). 



3 Numerical results 

PYQUEN partonic energy loss model was applied to cal- 
culate the nuclear modification factor of charged hadrons 
and the imbalance in dijet transverse energy in central 
PbPb collisions at y/s NN = 2.76 TeV for various scenarios 
of medium-induced radiation. The numerical results were 
compared with the LHC data applying the same kinematic 
cuts as in the experiments. 



(2) 3.1 Nuclear modification factor 



The nuclear modification factor Raa is defined as a ratio 
of particle yields in AA and pp collisions normalized on the 
number of binary nucleon-nucleon sub-collisions (N co n) : 



Raa(pt) 



d N I d"r\dp r p 
{N coU ) <PN^/dr)dp T 



(4) 



where px is the transverse momentum and r\ is the pseudo- 
rapidity. In the absence of nuclear effects (in initial or final 
states) for high px it should be Raa = 1- Figure Q] shows 
the calculated nuclear modification factor Raa for charged 
hadrons (|?7| < 0.8) as compared with ALICE data [5] in 
5% of most central PbPb collisions (T max =0.85 GeV for 
the left plot and 1 GeV for the right plot). Three options 
for the angular spectrum of in-medium gluon radiation 
(UJ, ^ and © were considered, collisional energy loss and 
nuclear shadowing being also taken into account. Higher 
initial temperatures result in more suppression of high- 
Pt hadrons, while a wider angular radiation results in a 
stronger dependence of Raa on px- The PYQUEN scenar- 
ios with T™ ax =1 GeV and wide- angular radiation (fJJ) or 
([3]) seem closer to the data than other options. However 
even in this case, the simulated p-r-dependence of -Raa 
looks weaker than the tendency of the ALICE data. The 
large systematic uncertainties of the data (due to the ab- 
sence of pp data at the same energy) still complicate a 
rigid quantitative interpretation of this observation. But 
if such discrepancy would be confirmed in the future with 
higher statistics and reference pp data at yfs — 2.76 TeV, 
in our opinion it could indicate that the real energy de- 
pendence of partonic energy loss AE is weaker than the 
expected one. It is worth to remind here that the energy 
dependence of BDMS radiative loss [2T1I22] (the dominant 
mechanism of medium-induced energy loss in PYQUEN) 
is between AE (x\aE and AE cx \f~E~ (depending on the 
energy of the radiated gluon and on the properties of the 
medium). 

Let us discuss some other model-dependent uncertain- 
ties of the results obtained. Excluding the collisional loss 
from the simulations increases Raa by ~ 15 -j- 20% for 
Pt — 10 4 20 GeV/c. Excluding the nuclear shadowing in- 
creases Raa by ~ 30 -j- 20% for the same kinematic range 
and weakens its px- dependence. The choice of the medium 
formation time tq within its reasonable range has rather 
moderate influence on the strength of medium-induced 
partonic energy loss. We have found that increasing tq 
by factor 2 (at fixed T max = 1 GeV) reduces Raa by 
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Fig. 1. The nuclear modification factor Raa for charged hadrons in 5% of most central PbPb collisions at \/s NN = 2.76 TeV. 
The points are ALICE data [6] (the error bars show the statistical uncertainties, the boxes show the systematic errors), the 
histograms are simulated PYQUEN events. The dotted, solid and dashed histograms correspond to the parametrizations of the 
radiation angular spectrum ([T}, ([2]) and Q respectively. The initial maximal temperature of quark-gluon matter T™ ax =0.85 
GeV (left) and 1 GeV (right). 
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Fig. 2. The event fraction of dijet events with the asymmetry ratio Aj (left) and with the leading jet energy (right) in 
10% of most central PbPb collisions at v^nn ~ 2.76 TeV. The points are CMS data [8] (the error bars show the statistical 
uncertainties), the histograms are simulated PYTHIA (dash-dotted) and PYQUEN (other styles) events. The dotted, solid and 
dashed histograms correspond to the parametrizations of the radiation angular spectrum UJ, ([2]| and Q respectively. The 
initial maximal temperature of quark-gluon matter T™ ax =l GeV. 



~ 13% without changing its pr-dependence. Note that the 
jet quenching effect gets stronger at larger To due to slower 
cooling of the matter, which implies a jet parton spending 
more time in the hottest regions of the medium. As a result 
the rescattering intensity increases (it starts little later for 
larger To, but this factor is much less significant). Fixing 
To = 0.2 fm/c and decreasing T™ 3 * to 0.85 GeV allow us to 
get the similar results for Raa(pt) as shown in Figure [T] 



(right). The jet quenching effect is also sensitive to the 
space-time evolution of a dense matter. For calculations 
presented, the medium was treated as a boost-invariant 
longitudinally expanding perfect quark-gluon fluid (well- 
known scaling solution obtained originally by Bjorken |31j , 
T(t) = ^(to/t) 1 / 3 ). In principle, other scenarios of QGM 
space-time evolution for PYQUEN are also envisaged [TH 
26 . For example, the presence of fluid viscosity slows down 
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the cooling rate, i.e., in fact the effective temperature 
of the medium gets higher as compared with the per- 
fect QGM case. On the other hand, the influence of the 
transverse expansion on the intensity of parton rescatter- 
ing is inessential for high initial temperatures T™ ax T c 
(T c ~ 200 MeV is the critical temperature). 



3.2 Dijet transverse energy imbalance 

Dijet asymmetry ratio Aj is defined [TIE] as a difference 
between transverse energies of two jets jl and j2 normal- 
ized on its sum: 

Aj = \ \ . (5) 

El 1 K ' 

In order to define calorimetric jets, we use the toy jet 
finding model. The jet energy was determined as the total 
transverse energy of the final particles collected around the 
direction of the leading particle inside the cone of radius 
R = ^Arj 2 + Aip 2 = 0.5 (77 and <p are the pseudorapidity 
and the azimuthal angle respectively). Then the final jet 
energy was smeared using Gaussian distribution with a 
proportional to the square root of a "true" jet energy. 
The same kinematic cuts as in the CMS analysis [8] were 
applied: |t7 j1 < 2 | < 2, \tp^ ~ ip> 2 \ > 2.1, E$ > 120 GeV, 

Ef > 50 GeV. Note that we compare our simulations with 
the CMS data having harder cuts on minimum jet energy 
than ATLAS data, in order to reduce complications facing 
jet reconstruction in high multiplicity background [32111 lj . 

Figure [2] (left) shows the calculated fraction of dijet 
events with the asymmetry ratio Aj for three options of 
the angular spectrum of in- medium gluon radiation ([T|), 
([2]) and ([3]) respectively. The choice of main parameters 
(T™ =1 GeV, t = 0.1 fm/c) was fixed from fitting the 
LHC data on Raa{pt) discussed in previous sub-section. 
The collisional energy loss is taken into account, the nu- 
clear shadowing effect being negligible for such high Et 
jets. PYQUEN results are compared with CMS data [8] 
for 10% of most central PbPb collisions and with PYTHIA 
simulations (in fact, PYQUEN without medium effects). 
The option with the small-angular radiative energy loss 
Q]) does not result in any additional dijet asymmetry 
as compared with PYTHIA (small additional asymmetry 
seen for dotted histogram is due to the collisional loss). 
Indeed, in this case the bulk of in-medium emitted gluons 
fly close to the parent parton direction and still belong 
to the jet. Therefore, the small-angular medium-induced 
radiation softens particle energy distributions inside the 
jet and increases the multiplicity of secondary particles, 
but does not affect the total jet energy [37] . On the other 
hand, the wide-angular radiation @ and © generates the 
significant dijet asymmetry due to the substantial frac- 
tion of "out-of-cone" jet energy loss. We would not ex- 
pect the perfect agreement of PYQUEN simulations with 
the data using our oversimplified treatment of jet finding 
(without taking into account realistic detector effects, jet 
reconstruction details and high multiplicity background 



fluctuations), but at least qualitative reproducing of the 
dijet asymmetry can be achieved with the wide-angular 
radiative and collisional partonic energy loss. This inter- 
pretation is in the agreement with the observation of CMS 
collaboration that the dijet momentum balance is recov- 
ered when integrating low transverse momentum particles 
distributed over a wide angular range relative to the di- 
rection of the away-side jet [5] . 

Figure [3] (right) represents the distribution over the 
leading jet energy EL . This distribution seems quite sta- 
ble. It does not depend significantly on jet quenching sce- 
nario, and is similar for PYQUEN, PYTHIA and CMS 
data. It confirms that the considered dijet configurations 
are dominated by a "surface emission" , when the dijet pro- 
duction vertex is close to the surface of the nuclear over- 
lapping area. In this case one partonic jet escapes from 
the dense zone almost without interactions and produces 
the leading hadronic jet, while a second partonic jet is af- 
fected by medium-induced multiple interactions and pro- 
duces the "quenched" hadronic jet. 

Note that the absolute number of dijet events (having 
leading jet with EL >120 GeV and away-side jet with 
>50 GeV) is suppressed in PYQUEN as compared 
with PYTHIA due to reducing the number of away-side 
jets observed above the threshold 50 GeV (it corresponds 
to increasing the "monojet" events). The dijet suppression 
factor depends on the radiation angular spectrum and is 
estimated as 0.53, 0.36 and 0.29 for PYQUEN options (JTJ), 
([2]) and ([3]) respectively. 



4 Conclusions 

The first LHC data on high-p^ hadron and dijet spectra in 
PbPb collisions at v^nn = 2.76 TeV are analyzed in the 
frameworks of PYQUEN jet quenching model. The pre- 
sented studies for the nuclear modification factor of high- 
Pt hadrons and the imbalance in dijet transverse energy 
support the supposition that the intensive wide-angular 
("out-of-cone") medium-induced partonic energy loss is 
seen in central PbPb collisions at the LHC. This inter- 
pretation is in the agreement with the recent observation 
of CMS collaboration that the dijet momentum balance 
is recovered when integrating low transverse momentum 
particles distributed over a wide angular range relative to 
the direction of the away-side jet. 

However, there are still some complications for a more 
quantitative interpretation of the jet quenching pattern at 
the LHC. The measured by ALICE nuclear modification 
factor i?AA is affected by large systematic uncertainties 
due to the absence of pp data at the same energy. More- 
over, the contribution of nuclear shadowing to Raa seems 
visible at least up to pr ~ 25 GeV/c and should be taken 
into account. Concerning ATLAS and CMS data on the 
dijet asymmetry, simple generator level studies are likely 
not enough to quantify this observable accurately (jet find- 
ing details and high multiplicity background fluctuations 
may be important). 
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